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THE  GENERATION  AND  RECOGNITION  OF  FLICKER  NOISE 


James  A.  Barnes 
Abstract 

Flicker  noise  is  defined  as  noise  whose  power  spectral 
density  varies  inversely  proportional  to  the  Fourier 
frequency.  Flicker  noise  appears  to  be  present  in  every 
electronic  devicq  but  comprehensive  explanations  as  to  its 
source  and  level  are  lacking.  It  is  shown  that  one  can 
approximate  flicker  noise  by  analog  methods  which  are 
amenable  to  the  use  of  analog  computers.  It  is  also  shown  that 
certain  diffusion  processes  and  certain  Poisson  processes  can 
generate  flicker  noise.  A method  of  generating  flicker  noise 
using  a digital  computer  is  presented. 

With  a little  practice  one  can  learn  to  recognize  noise 
processes  with  various  power  law  spectra--some  examples 
are  shown.  A review  of  techniques  for  inferring  the  power 
spectral  density  from  time  domain  measurements  which  have 
been  developed  recently  by  several  authors  is  presented. 


Key  W ords: 


Noise,  Flicker  Noise,  Diffusion,  Poisson  Process, 
Spectral  Density,  i-  noise,  Effective  Stationarity. 
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Introduction 


Flicker  noise  is  the  name  given  to  a random  noise  whose  power 
spectral  density  varies  as  over  a very  large  spectral  region. 

Equivalently  in  the  time  domain,  a time  series  whose  standard  deviation 
is  independent  of  the  sampling  time  is  a flicker  noise.  While  it  is  possible 
to  present  compelling  arguments  why  this  behavior  cannot  persist  to  zero 
Fourier  frequency,  experiments  to  determine  a low  frequency  break  in 
behavior  have  failed  for  most  (if  not  all)  devices.  Measurements  have 
indicated  a continuation  of  this  behavior  at  one  cycle  per  year  in  some 
devices  [Atkinson,  et  al.  , 1963]. 

Part  of  the  difficulty  in  analysis  of  this  type  of  noise  may  be  ascribed 
to  its  very  mild  divergence  characteristic.  That  is,  if  one  assumes  that 
an  oscillator's  frequency  is  modulated  by  a flicker  noise,  and  one  can 
measure  the  frequency  with  an  uncertainty  of,  say,  ±A  in  five  milliseconds 
then  one  can  show  that  the  oscillator  will  be  within  about  ±5A  of  that  fre- 
quency at  an  age  of  the  Universe  later!  This  is  a logarithmic  divergence 
in  the  time  domain,  and  it  is  so  mild  as  to  give  experimenters  a false 
sense  of  security  and  to  cause  them  to  consider  the  process  to  be  stationar 
It  is  interesting  to  note  that  the  false  assumption  of  stationarity  can,  in 
specific  situations,  yield  correct  answers  (see  Appendix). 

The  present  paper  places  emphasis  on  noise  processes  which  have  a 
power  spectral  density  varying  as  |w|-1.  However,  it  is  of  value  to 
consider  this  noise  in  its  relation  to  noises  whose  power  spectral  densities 
are  proportional  to  | oo  | for  constant  values  of  a. 

Implicit  in  most  of  the  discussions  in  this  paper  is  the  assumption 
that  one  may  meaningfully  consider  power  spectral  densities  of  the  form 

= h |co|  • 

This  is  not  a trivial  assumption  because  for  a 2-  1 one  can  show  that  <f>(t) 
is  not  covariant  stationary  [Barnes  and  Allan,  1966].  Thus  the  assumption 


that  S 

There  are  given  in  this  paper  two  plausibility  arguments  to  support  the 
assumption  that 


- a 


(to)  = h [to  I raises  doubt  that 


can  even  be  defined  for  cf)(t). 


S^co)  = h 


O'  > 1 


is  a meaningful  relation. 

First  and  foremost,  experimentally  one  encounters  noise  processes 
which  appear,  in  fact,  to  be  described  by  such  power  spectra  and  very 
often  o'  S:  1 . Here  one  may  well  wish  to  compare  the  actual  experimental 
methods  with  theory. 

Secondly,  it  is  entirely  possible  that  these  power  spectra  break  into 
a well  behaved,  convergent  spectral  law  as  to  — 0.  Such  cutoffs  are  at 
extremely  low  frequencies  for  some  devices  since,  in  general,  these 
cutoff  frequencies  have  not  been  observed.  One  may  thus  assume  a lower 
cutoff  frequency,  a>^,  and  construct  the  theory  accordingly.  Since  has, 
to  date,  not  been  an  observable  in  most  cases,  one  considers  only  those 
quantities  which  are  insensitive  to  for  very  small  co^.  This  could  be 
called  a cutoff  independent  formulation.  This  method  of  formulating  the 
problems  of  power  spectra  with  the  specific  form 


S^co)  ^ h jojjf  + jco  | 

is  considered  in  some  detail  in  the  Appendix.  One  can  show  [Cutler  and 
Searle,  1966]  that  many  quantities  of  interest  are  not  critically  dependent 
on  the  form  of  for  u>  smaller  than  the  reciprocal  observation  time-- 

indeed,  this  is  precisely  why  co^  is  not  normally  an  observable. 
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1.  The  Generation  of  Flicker  Noise 
1.  1 Analog  methods 

An  operational  amplifier  is  a device  whose  transfer  function,  g(co), 
is  approximated  by  a large,  real,  negative  number,  -K,  over  a very  large 
spectral  region  (0  ^ w ^ gj-^)  which  is  assumed  to  include  the  frequency 
domain  of  interest.  If  one  connects  such  an  operational  amplifier  with  the 
two  impedances,  and  Z^,  as  shown, 


the  overall  transfer  function,  G(oo),  is  given  by 


G(w)  = 


eo<“> 

e,  (u) 


-K 


1 + (1  + K)- 


Thus,  if  K » 


Consider  the  case  where  the  first  impedance  is  a resistor,  Z^  = R, 

and  the  second  a "fractional  capacitor"  [Carlson  and  Halijak,  1964], 

2,  = j — - — | > where  Z is  defined  as  the  impedance  of  a — - th  "fraction" 

2 \jcoc  / 2 2 

of  a capacitor.  Actual  networks  which  can  realizre  impedances  of  this  form 
are  considered  below.  The  transfer  function  of  this  device  then  becomes 
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For  convenience  and  simplicity  it  is  desirable  to  follow  this  device 
by  an  inverting  amplifier  to  give  an  overall  transfer  function  of 


g'M 


(0  K 


(i-i) 


If  the  input  signal,  e^,  has  a power  spectral  density,  S.(u>),  and, 
similarly,  the  output  signal,  e^,  has  a power  spectral  density,  Sq(oo), 
then 


S («)  = |g'(u>)  1 2 S.  (co). 

o 1 


(1-2) 


For  a white  input  noise  of  power  spectral  density  h and  t = 1 second, 

Sq(cu)  ^ hlcol"0,  (1-3) 

for  M^K2  > > 1.  Thus,  one  may  generate  | co  | ^ noise  from  a white 
noise  source,  an  operational  amplifier,  and  a "fractional  capacitor.  " 

It  is  of  interest  to  determine  the  operation  of  this  device  in  the  time 
domain,  also.  The  impulse  response  function  of  a filter  is  the  Fourier 
transform  of  the  transfer  function.  For  a filter  whose  transfer  function 


4 


a 


Foster, 


one  can  show  that  the  impulse  response  function  is  [Campbell  and 
1948  (Pair  516)] 


- 1 


Thus,  if  e.(t)  is  the  input  signal  in  the  time  domain,  the  output,  e (t),  is 
1 o 

the  convolution  of  e.(t)  with  the  impulse  response  function.  That  is. 


eQ(t) 


e.(t  )dt  , 
l o o 


(1-4) 


or  e (t)  is  the  - fold  integral  [Courant,  1957]  of  e.(t).  Note  that  for 

(X 

“ = 1 , one  has  a 1 ’whole1 1 capacitor  and  the  output  is 


eo(t) 


e.(t  )dt 
1 o 


o' 


(1-5) 


This  is  the  proper  result  for  an  integrator  as  used  in  analog  computers. 

Typically,  integrators  used  in  analog  computers  realize  a K ^ 10  . Thus 

one,  theoretically,  has  an  integrator  working  down  to  ~ one  cycle  per  12 

days.  There  are,  of  course,  other  effects  which  limit  its  operation  long 

before  this  low  frequency  is  actually  realized. 

a 

~ 2* 

1.  Z Realization  of  the  impedance  (jco) 


O.  Heaviside  [1950]  noted  that  the  impedance  of  a semi-infinite  lossy 


transmission  line  has  the  form 


T jgj 


One  may  show  that  this  is  reason- 


able by  considering  the  lumped-constant  transmission  line  shown. 
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p6X 

AAA 


P6X 

aaA 


P6X 

AAA- 


C6X 


c6X 


Since  the  line  is  infinite  in  extent  and  has  impedance  Z,  nothing  is  changed 
by  adding  one  more  section* 


p6X 


The  solution  of  this  network  is 


Z2  = -A-  + Zp6X. 

JGJC 


In  the  limit  as  6X 


(a  continuous  line). 


(1  -6)  reduces  to 


(1-6) 


7 2 — 

jojc 


where  p and  c are  the  resistance  and  capacitance  per  unit  length  of  the 
line. 

— ct 

Additional  methods  of  constructing  impedances  of  the  form  (jcj) 
may  be  found  in  the  literature  [Carlson  and  Halijak,  1961;  Lerner,  1963; 
and  Roy  and  Shenoi,  1966].  An  intuitive  and  relatively  simple  method 


6 


of  constructing  such  impedances  was  devised  by  the  author.  Consider 
the  network: 


The  transfer  function  of  this  network  is 

1 + jtOT 


gM  = 


1 + ja>(Tl  + T2) 


(1-7) 


where  t ^ = R^C  and  t ^ = R^C.  For  smaH  frequencies  (u>(t  i + T 2)  < < *)  » 

g(a>)«sl;  for  large  frequencies  (cot  > > l)  , g(c 0)  ^ 7 . 

2 Tl+T2 

If  one  has  a sequence  of  such  filters  such  that  for  the  i-th  filter 

(i)  _ (o) 

T l = (P)  T 1 

and 

t ^ = (PlK  2°)  for  i = 1, 2,  • • • N,  and  (3  < 1, 

and  these  filters  are  cascaded  with  appropriate  isolation  amplifiers 
between  successive  filters,  the  over-all  transfer  function  becomes 


G(o>)  = 


N 


i=l 


l + jcOT^fp)1 


1 + 


+ T 2°^) 


(1-8) 


At  frequencies  near  where  got  ^°^(p)^sa  1>  the  transfer  function  has 


approximately  the  value 


G (us)xd 


(o) 

1 


(o)  (O) 

T1  +r2 


N 
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At  frequencies  near  ojt  1 the  transfer  function  has  approximately 

the  value  unity.  Thus,  over  an  interval  of  frequency  from 

oo , = - to  oo„  = . r the  (voltage)  transfer  function  is  reduced  by 


1 (o) 

T1 

the  factor 


N (o) 


P T 


1 


(O) 

1 


N 


(O)  x (o) 
‘T1  +T2 


If  this  filter  is  to  approximate  a filter  whose  transfer  function 
(in  magnitude)  varies  according  to  the  relation 


|g'(u) 


O' 
, 2 


00 


then 


For  flicker  noise,  a = 1 and  (1-9)  becomes 


(1-9) 
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Setting  N = 4,  and  t |°^  + r ^ ^ = 3t  |°\  one  has  (3  = and  (1  -8)  becomes 


/ x _ /_S_j_9_\  /S+  81  \ / 

G(co)  “ V3S+  1/  \3S  + 9/  V3S+  81/ V 


S+  729 
3S+  81/  \3S  + 7Z9 


(1-10) 


where  S = ju>T 


(o) 


am 


Equation  (1-10)  is  the  transfer  function  of  an  active  filter  (with  isolation 
plifiers  between  sections)  which  should  approximate  2 over  a 


relative  frequency  range  of  ~ 9 or  nearly  4 decades  of  frequency.  If 

/ A \ ~ 

one  were  to  consider  (1-10)  as  a function  which  approximates  (-t-j-J2,  one 
could  consider  an  impedance,  Z(oo),  which  has  exactly  this  same 
functional  form.  That  is. 


/S+i  \ 

(S+9\ 

/S+  81  \ / S+  7Z9  \ 

\3S+  1/ 

\3S  + 9/ 

V3S+  81/  V3S+  7Z9/ 

(1-11) 


One  can  decompose  (1-11)  into  a realizable  impedance: 
Z.  Z50  0.  96Z  0.  3Z3  0.  0863 


0. 0Z34 


By  scaling  the  resistor  and  capacitor  values,  one  obtains  a circuit 

nr 

whose  impedance  approximates  y-. — for  A = 1816  , within  ±0.  5 dB 

from  about  1 Hz  to  4 kHz.  By  using  such  an  impedance  with  an  operational 

amplifier  as  in  Section  1.  1 (impedance  Z ),one  obtains  an  amplifier  whose 

' A 4 Z 


(magnitude)  gain  differs  from 


as  shown  in  Figure  1. 
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CQ  CQ 

xi  o T3 
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The  circuit  diagram  of  the  impedance  used  for  Figure  1 is  shown  below. 


4.  475  1.  913  0.  64Z4  0.  1716 


Also  of  interest,  if  one  intercharges  Z and  Z using  the  previous 

1 z +1 
filter,  one  obtains  a filter  with  a transfer  function  proportional  (joo)  2. 

Thus,  if  one  applies  a flicker  noise  signal  to  the  input,  the  output  is 

white  noise. 

1.  3 Flicker  noise  in  Diffusion  Processes 

Heaviside's  semi-infinite  lossy  transmission  line  may  be  considered 
the  electrical  analog  of  a diffusion  process.  Consider  a semi-infinite  rod 
which  obeys  the  diffusion  equation 


Ziiijkl)  . a2  MGk!>  _ . 4 Tr  p(X,  t)  (1-12) 

5X2  St 

where  ip(X,  t)  might  be  temperature  or  an  impurity  concentration, and 
p(X,  t)  is  the  source  function,  e.  g.  , a heat  source.  The  Green’s  function 
[Morse  and  Feshbach,  1953]  for  this  problem  is 


gfc-Xo.i-t,)  = 


e 


a2  {X"Xo)2 
4(t-t0) 


u(t-to) 


where  U (t  ) is  the  unit  step  function  defined  by  the  relation 


U(t) 


1>  T ^0 

0»  T < 0. 
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Thus,  if  the  source  for  this  system  is  e.(t)  located  at  X = 0,  the  end  of  the 
rod,  the  solution  of  (1-12)  may  be  written 


(1-13) 


At  the  end  of  the  rod,  X = 0,  (1-13)  is  formally  identical  to  (1-4)  for  a = 1 


lp[0,t)  is  a flicker  noise  process. 

An  obvious  implication  is  that  impurity  diffusion  in  a semiconductor 
should  cause  a flicker  noise  process  at  the  surface--a  critical  region  for 
semiconductor  devices.  It  is  interesting  to  note  that  flicker  noise  is  often 
referred  to  as  "semiconductor  noise.  " The  apparent  lack  of  a strong 
temperature  dependence  of  actual  flicker  processes,  however,  does  not 
seem  to  support  the  present  model  of  impurity  diffusion  as  the  cause  of 
flicker  noise.  A different  diffusion  model  has  been  considered  by 
J.  M.  Richardson  [ 1 9 5 0 ] and  shown  to  give  flicker  noise  and  yet  not  be 
significantly  temperature  sensitive.  However,  some  assumptions  in  the 
Richardson  model  do  not  seem  entirely  physical. 

1.  4 Quantum  Mechanical  Model 

The  treatment  here,  due  primarily  to  Dr.  D.  Halford  [1967]  , 
is  based  on  discrete  events  in  a Poisson  process  and  hence  is  easily 
amenable  to  quantum  statistics. 

Consider  a real,  "localizable"  function  h(x).  For  convenience,  let 


(except  for  the  factor  - — ),  and  thus,  if  e,(t)  is  a white  noise  process, 

H.  1 


(1-14) 
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It  is,  thus,  true  that  h(X)  has  a Fourier  transform,  g(y),  defined  by  the 
relations 


and 


(1-15) 


h(X)e'lXydX. 


(1-16) 


One  requires  h(x)  to  be  a "reasonable"  function  in  the  sense  of  being  nearly 
realizable  in  some  physical  device.  Therefore,  one  further  assumes  h(X) 
is  finite  for  all  values  of  X. 

It  is  true  that,  for  y s 0,  (1-16)  may  be  written  in  the  form 


|g(0)|=|  f h(X)dX  | 

•Coo 


< 


(1-17) 


by  virtue  of  (1-14).  Thus,  in  the  region  of  small  y,  if  one  were  to 
approximate  the  energy  spectrum. 


U(y)  $ 1 g (y)  1 2» 


by  a power  law. 


one  would  have 


U (y)  =» 


(1-18) 


for  "small"y,  and,  hence, 

Y £ 0. 

o 


(1-19) 
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In  the  region  of  "large"  y,  the  total  energy  must  also  be  finite;  i.  e.  , 


g(y)  |2dy  < °°,  (l -20) 

where  y is  a "large"  number.  This  implies  that  if  U (y)  asymptotically 

1 - ,Y 

approaches  the  function  B |y  | “ as  y -*  00 , 

y < - i. 

co 

It  can  be  shown  that  the  conditions  on  h(x)  are  sufficient  to  guarantee  that 

y z - 2.  (1-21) 

00 

While  it  is  possible  to  let  b(X)  be  a very  general  function,  it  is 
sufficient  for  the  present  treatment  to  further  restrict  the  acceptable 
class  of  functions,  h(X).  In  particular,  define  y (y)  such  that 

U (y)  = A | y | Y for  all  y. 

Only  those  functions,  h(X),  will  be  considered  for  which  "Y  (y)  is  a 
monotonic  function  of  y from  y to  Y as  y goes  from  0 to  00  . This  avoids 

O oo  ' ° 

some  problems  in  defining  effective  times  of  the  functions  generated  by 
h(X). 

One  may  now  consider  a Poisson  process  composed  of  a superposition 
of  functions  obtained  from  h(X)  by  scaling  amplitude  and  time  and  also  by 
translating  time;  i.  e.  , 

m.  (t)  = a.h 
i i 


jr) 


(1-22) 
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Since  this  is  assumed  to  be  a Poisson  process,  the  distribution  of 
occurrence  times,  t.,  is  random,  uncorrelated  and  uniformly  distributed 
over  the  time  interval  of  interest. 

Let  the  t.  's  and  a.'s  be  distributed  according  to  some  probability 
distribution  density  function,  p(T,a).  That  is,  the  probability  of 
observing  a "life-time"  between  t and  t + dr  and  an  amplitude  between 
a and  a + da  is 

p(r  , a)dT  da. 


The  Fourier  transform  of  iru(t)  may  easily  be  shown  to  be 


Vi (W)  V mi 


(t)e-lwtdt 


(1-23) 


-icot. 

= aiT  ie  g (cot  i) , 

where  g(y)  is  defined  in  (1-15)  and  (1-16).  The  energy  spectrum 
associated  with  m.(t)  is  then 

l 

b.Hl2  = aiTi  lg(“T.)|2. 

If  R is  the  average  rate  of  occurrence  of  events,  m.(t),  then  the 
average  rate  of  energy  expenditure  (i.  e.  , power)  at  frequency  go  due  to 
events  of  life-times  between  t and  t + dT  and  of  amplitudes  between  a 
and  a + da  is  given  by 


d2  S (go)  = Ra2T  2p  ( t , a)  | g (got  ) | 2 dadT  . 


Thus  the  total  power  spectrum  due  to  all  events  is 


CD  _ 00 


>// 


S(go)  = R / / a2  T 2 p (t  , a)  | g (coT  ) 1 2dadT  . 


(1-24) 
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Consider  the  function  A2(t)P(t)  defined  by  the  relation 


A2  (t)P(t) 


CD 

a2p(T,  a)da} 


and  P(t)  is  the  probability  density  for  t without  regard  to  amplitude  and 
A2  (t)  is  the  mean  square  amplitude  for  the  actual  events  of  lifetime  t. 
Equation  (1-Z4)  can  now  be  written  in  the  form 


S(co)  = R 


t2A2 (t)P(t)  |g(u-r)2  I dr. 


d-25) 


If  it  should  happen  that 

t2A2(t)P(t)  = Bt^  (1-26) 

over  some  range  of  t (say  T < t<  t^)  and  be  zero  (or  practically  so) 
outside  this  range,  then  (1-25)  becomes 


S(u)  = 


|g(y)  l2y6dy- 


(1-27) 


If  it  is  true  that  y = cot  and  y = cot  reasonably  include  the  major 

O O « oo 

portion  of  the  area  under  the  function  |g(y)  |2y6  for  positive  y,  then  (1-27) 
may  be  approximated  by 


s(w)  « | *jf+  6 J |g(y)l2y6dy 

0 


and  for  6 = 0 one  obtains 


S(co)  & 


RBtt 


CO 


(1-28) 


16 


by  virtue  of  (1  - 1 4),  (1  - 1 5),  and  (1-16)  and  the  symmetry  of  | g (y)  |2.  Within 
the  restrictions  noted  above,  (1-28)  is  valid  and  hence  the  process  generates 
a flicker  noise.  These  restrictions  on  a)  may  be  expressed  in  the  following 
forms  for  6 = 0: 


1 


Ci)  T 


g(y)  I2dy<<  i 


(i)  T. 


(1-29) 


or 


I g(y)  l2dy  < < i. 


00 


I g (y)  1 2 dy  < < 


1. 


(1-30) 


Obviously  these  conditions  cannot  be  satisfied  at  gj  = 0. 

More  extensive  and  complete  solutions  of  the  integral  equation 


h |a> 


-a 


2A2(t)P(t)  |g(wr)  |2dT 


have  not  been  considered  in  the  literature  to  this  author's  knowledge. 

1.  5 Comments  on  Current  Papers 

For  6 = 0,  (1-26)  may  be  written  in  the  form 

A2(t)P(t)  =JL.  (1-31) 

2 
T 

That  is,  a necessary  condition  for  the  generation  of  flicker  noise  by  the 
process  just  described  is  that  (1-31)  is  valid  over  a sufficient  region. 
Unfortunately,  one  can  find  in  the  literature  [van  der  Ziel,  1950; 
Mueller,  1965;  and  Jordan  and  Jordan,  1965]  this  condition  reduced 
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nearest  integer.  Starting  at  location  t.  , the  integer  1 is  added  to  every 
location  up  to  and  including  a location  congruent  to  (t.  + T.)  mod  T.  By 
using  modular  arithmetic  in  the  storage  process  one  insures  that  the 
entire  data  are  equivalent  to  a segment  of  an  infinite  periodic  process. 

An  alternative  is  to  limit  t.  + t\  ^ T and  use  only  the  latter  portion  of 
the  data. 

Repeating  the  above  process  a total  of  4096  times  with  T = 10Z4, 
yielded  the  results  shown  in  figure  2.  Also  included  in  figure  2 are  the 
(finite)  integrals  of  the  flicker  noise  samples.  Using  methods  described 
in  Section  2,  six  independent  samples  of  flicker  noise  generated  by  this 
scheme  were  analyzed.  The  data  were  consistent  with  a power  spectral 
density  of  the  form 

S(oo)  = h | oo  | 01 

where  a was  determined  to  be  0.  98±.  02.  The  random  numbers  used  for 
these  tests  were  computer  generated  by  a pseudorandom  process  [Andrew, 
1966;  and  Chambers,  1967]. 
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Figure  2 
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2.  Noise  Recognition 


2.  1 General  Remarks 

The  most  comprehensive  treatment  of  recognition  of  noises  with 
the  power  law  types  of  power  spectral  densities  is  to  be  found  in  the 
Proceedings  of  the  IEEE,  Vol.  54,  No.  2,  February  1966.  Of  particular 
note  are  the  papers  by  Cutler  and  Searle  [ 1 966],  Vessot,  et  al.  [ 1 96 6 ] , 
Allan  [1966],  and  Barnes  [1966].  These  papers  deal  with  the  inference 
of  power  spectral  densities  from  certain  time  domain  analyses.  While 
one  can  approach  the  classification  problem  directly  in  the  frequency 
domain  either  by  analog  techniques  or  other  analysis  [Blackman  and 
Tukey,  1958],  one  is  normally  interested  in  predicting  performance  both 
in  the  frequency  domain  and  the  time  domain.  The  time  (frequency) 
domain  implications  of  the  frequency  (time)  domain  behavior  are  of  great 
value. 

2.  2 Visual  Inspection 

One  of  the  most  surprising  aspects  of  noise  characterized  by  a 
power  law  type  of  power  spectral  density  is  the  fact  that,  with  a little 
practice,  a visual  inspection  of  the  data  is  often  adequate  to  classify  the 
spectral  type.  Of  course,  one  often  has  other  information  about  the  noise 
source  and  hence  often  one  has  reason  to  suspect  one  type  of  noise  over 
another  before  the  data  are  acquired. 

It  is  apparent  that,  if  one  has  a noise  whose  spectral  density  is  of 
the  form 

S (u ) = h | co  | a, 

then  for  a ^ 1 the  process  is  high  frequency  divergent.  One  certainly 
does  not,  in  reality,  observe  such  a process  since  the  measuring 
apparatus  has  a finite  high  frequency  limit.  What  one  does  observe, 
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however,  is  that  for  processes  where  a < 2,  the  visual  appearance  can  be 
drastically  changed  by  the  high  frequency  cutoff. 

While  it  is  true  that  for  a ^ 1 the  process  is  low  frequency  divergent, 
Cutler  [1966]  points  out  that  a finite  data  sample  analysis  is  equivalent  to 
high-pass  filtration.  Thus,  the  low  frequency  divergent  character  is 
obscured  to  some  extent.  Indeed,  the  operation  of  "zeroing"  amplifiers, 
oscilliscopes,  and  various  measuring  devices,  is,  to  some  extent,  high- 
pass  filtration. 

Thus  for  valid  noise  comparisons,  one  should  compare  those  noise 
processes  with  similar  products  of  data  length  with  high  frequency  cutoff. 

In  figure  3 a number  of  noise  samples  which  were  computer  generated  are 
shown.  The  cutoff  (cycle)  frequency  data  length  product  is  about  500  for 
figure  3. 

2.  3 Standard  Deviation  Measurements 

A familiar  statistic  of  a noise  process  is  the  standard  deviation  of 

different  values  of  the  noise,  X(t).  Normally  when  considering  a continuous 

noise  process  one  averages  the  noise  value  over  some  time  interval,  t, 

to  obtain  a finite  set  of,  say,  N measurements.  As  pointed  out  above, 

for  some  noises  the  system  bandwidth,  co^,  is  a significant  parameter. 

If  there  exists  "dead  time",  t^,  between  successive  averages  of  X,  then 

this  dead  time  may  also  be  of  significance.  Thus  one  may  compute  a 

quantity  s2(N,  t,  t ,, co , ) defined  by  typical  equations  for  the  sample  standard 
d h 

deviation: 

s=(N,T,Td.V  =— -Ztf.  (2-1) 

where  is  the  average  of  X(t)  over  the  n-th  time  interval  of  duration  t . 

Normally  one  is  interested  in  a sample  standard  deviation  as  an 
estimate  of  the  "true"  standard  deviation  which  would  be  obtained  in  the 
limit  as  N approaches  infinity.  For  many  of  the  types  of  noises  being 
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considered  here,  however,  this  limit  does  not  exist  (i.  e.  s diverges  as 
N -*  °°).  For  the  work  here,  the  "true"  standard  deviation  m N,  t,  t gj  ) 

V d h/ 

is  defined  as  the  limit 


M-  02 


M 

ct2(n,  t, Td>“h) = 5„im_  (' 

i=  1 


— V sf  (N,  t,  TJf 


where  s.  ^N,  t,  t^,  is  the  sample  standard  deviation  of  the  i-th  inde- 
pendent set  of  measurements  characterized  by  the  parameters  N,  t,t,, 

d 

and  CD  . This  approach  extends  the  utility  of  the  concept  of  the  standard 
h 

deviation  to  noises  whose  power  spectral  densities  may  be  described  as 


Sx(u)  = h|w 


- a 


where  a < 3.  For  noise  processes  where  a 2.  3 other  techniques  must  be 
used.  Since  one  can  never  actually  pass  to  the  limit  M -»  00 , one  measures 
again  a sample  standard  deviation. 

It  has  been  shown  by  Vessot,  et  al.  [1966]  and  Allan  [1966]  that  if 


Sx(w)  = h | to 


- a 


then 


ffa(N>T>Td>“h)  =kTM 


where  too  >>1,  kisa  function  of  N,  r u , and  \i  and  o-  are  related  as 
h d h 

shown  in  figure  4. 

It  is  thus  possible  to  plot  In  ct2(n,  t , t^,  co^  vs.  In  t and  infer  a from 
the  slope  of  this  plot  and  figure  4. 
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Figure  4 /j  - a Mapping 


2& 


For  co  - 1,  Allan  [ 1 966]  shows  that 


(j2(N,T,Td,a>h)=  A 


N-l  Jl+2 

l x (N  - n)  (nr)r 


N(N  - 1) 


n=  1 


t + T 


where  r = 


t to  > > l,and  A depends  only  on  to  , the  noise  level, 
h h 


and  noise  spectral  type. 


2.  4 The  Chi  Test 

Following  Allan  [1966],  one  may  define  the  quantity,  X,  by  the  relation 


X(N.a-)  = 


CT2(N,T,0,wh) 


on  2,  r,  o,  to 


(2-2) 


This  should  not  be  confused  with  the  statistician's  chi- squared  test.  . 

Provided  t<o  >>  1,  X(N,o')  is  a function  only  of  N and  the  spectral  type 
h 

of  noise  being  considered.  Table  1 lists  values  of  X(N,  a)  for  typical 
values  of  N and  a. 

2.  5 Finite  Differences 

The  first  finite  difference  of  the  quantity  X(t)  may  be  defined  as 

M.  = X(t  +T)  - X(t.). 

Similarly,  the  second,  third,  and  fourth  differences  are 


A2X.  = X(t  +2t)  - 2X(t  +T)  + X(t  ) 

A3X.  = X(t.  + 3t)  - 3X(t.  + 2t)  + 3X(t.  + t)  - X(t.) 
i i i l i 

and 

A4X.  = X(t.  + 4t)  - 4X(t.+  3t)  + 6X(t  + 2t)  - 4(t  + t)  + X(t.). 
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N 


tk 

4 

8 

16 

32 

-2.  0 

. 833 

. 750 

. 708 

. 687 

-1.  9 

. 845 

. 765 

. 724 

. 704 

-1.  8 

. 858 

. 782 

. 743 

. 722 

-1.  7 

. 871 

. 801 

. 763 

. 743 

-1.  6 

. 886 

. 822 

. 786 

. 767 

-1.  5 

. 902 

. 844 

. 812 

. 793 

-1.  4 

. 919 

. 870 

. 841 

. 824 

-1.  3 

. 937 

. 897 

. 873 

. 859 

-1.  2 

. 956 

. 928 

. 910 

. 899 

-1.  1 

. 977 

. 962 

. 952 

. 946 

-1.  0 

1.  000 

1.  000 

1.  000 

1.  000 

- . 9 

1.  023 

1.  041 

1.  054 

1.  062 

- . 8 

1.  049 

1.  088 

1.  116 

1.  136 

- . 7 

1.  077 

1.  139 

1.  188 

1.  223 

- . 6 

1.  106 

1.  197 

1.  270 

1.  327 

- . 5 

1.  138 

1.  261 

1.  365 

1.  450 

- . 4 

1.  171 

1.  332 

1.  475 

1.  598 

- . 3 

1.  208 

1. 412 

1.  604 

1.  777 

- . 2 

1.  247 

1.  501 

1.  753 

1.  993 

- . 1 

1.  288 

1.  602 

1.  928 

2.  257 

. 0 

1.  333 

1.  714 

2.  133 

2.  580 

. 1 

1.  381 

1.  840 

2.  374 

2.  978 

. 2 

1.  432 

1.  981 

2.  658 

3.  470 

. 3 

1.  487 

2.  141 

2.  993 

4.  082 

. 4 

1.  546 

2.  320 

3.  390 

4.  846 

. 5 

1.  609 

2.  522 

3.  862 

5.  802 

. 6 

1.  677 

2.  750 

4.  424 

7.  005 

. 7 

1.  749 

3.  007 

5.  093 

8.  523 

. 8 

1.  827 

3.  298 

5.  893 

10.  446 

. 9 

1.  910 

3.  627 

6.  851 

12.  888 

1.  0 

1.  999 

4.  000 

8.  000 

16.  000 

1.  1 

2.  095 

4.  421 

9.  379 

19.  974 

1.  2 

2.  198 

4.  900 

11.  040 

25.  062 

1.  3 

2.  308 

5.  442 

13.  041 

31.  593 

1.  4 

2.  426 

6.  059 

15.  457 

39.  992 

1.  5 

2.  552 

6.  759 

18.  376 

50.  815 

1.  6 

2.  687 

7.  554 

21.  908 

64.  788 

1.  7 

2.  832 

8.  460 

26.  187 

82.  855 

1.  8 

2.  988 

9.  490 

31.  377 

106.  253 

1.  9 

3.  154 

10.  663 

37.  677 

136.  597 

2.  0 

3.  333 

11. 999 

45.  333 

176.  000 

TABLE  OF  X VALUES 
TABLE  1 
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N 


64 

128 

256 

512 

1024 

. 677 

. 671 

. 669 

. 667 

. 667 

. 693 

. 688 

. 685 

. 684 

. 683 

. 712 

. 706 

. 704 

. 702 

. 702 

. 733 

. 727 

. 725 

. 723 

. 723 

. 756 

. 751 

. 748 

. 747 

. 746 

. 784 

. 779 

. 776 

. 774 

. 774 

. 815 

. 810 

. 807 

. 806 

. 805 

. 851 

. 847 

. 844 

. 843 

. 842 

. 893 

. 889 

. 887 

. 886 

. 886 

. 942 

. 940 

. 938 

. 938 

. 937 

1.  000 

1.  000 

1.  000 

1.  000 

1.  000 

1.  068 

1.  072 

1.  074 

1.  075 

1.  076 

1.  150 

1.  159 

1.  165 

1.  168 

1.  171 

1.  249 

1.  266 

1.  278 

1.  286 

1.  291 

1.  369 

1.  400 

1.  422 

1.  437 

1.  447 

1.  517 

1.  568 

1.  606 

1.  634 

1.  655 

1.  700 

1.  782 

1.  847 

1.  898 

1.  937 

1.  928 

2.  058 

2.  167 

2.  257 

2.  332 

2.  215 

2.  417 

2.  598 

2.  758 

2.  899 

2.  580 

2.  892 

3.  190 

3.  472 

3.  736 

3.  047 

3.  527 

4.  015 

4.  508 

5.  004 

3.  649 

4.  384 

5.  183 

6.  045 

6.  973 

4.  431 

5.  554 

6.  857 

8.  362 

10.  097 

5.  454 

7.  166 

9.  290 

11.  916 

15.  156 

6.  800 

9.  407 

12.  866 

17. 444 

23.  496 

8.  583 

12.  547 

18.  177 

26.  157 

37.  456 

10.  957 

16.  982 

26.  141 

40. 046 

61.  139 

14.  135 

23.  285 

38.  181 

62.  405 

101.  779 

18.  407 

32.  301 

56.  521 

98.  717 

172.  209 

24.  177 

45.  265 

84.  639 

158.  134 

295.  300 

32.  000 

64.  000 

128.  000 

256.  000 

511.  999 

42.  643 

91.  192 

195.  211 

418.  126 

895.  897 

57.  172 

130.  818 

299.  864 

688.  060 

1579.  703 

77.  064 

188.  766 

463.  475 

1139.  479 

2803.  483 

104.  373 

273.  778 

720.  174 

1897.  334 

5002. 701 

141. 955 

398.  853 

1124.  206 

3174.  015 

8969.  196 

193.  791 

583.  344 

1761.  937 

5331. 243 

16146.  Oil 

265.  432 

856.  107 

2771.  095 

8986.  198 

29168.  198 

364.  630 

1260. 222 

4371.  641 

15193.  658 

52856. 122 

502.  214 

1860.  075 

6915.  323 

25758.  904 

96042.  085 

693.  333 

2752.  000 

10965.  336 

43776.  034 

174933.  41 

TABLE  OF  X VALUES 
Table  1 (continued) 


x=^^;<03,N>T)>=aM  NjN^l)  |t|M 

\CT  (Z>T)/ 


. - a 


and  S(w)  = h |u>  | ; fi  - a - 1 , - 2<  fj,  < 2 
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For  Fourier  frequencies,  co,  such  that 


ut  > > 1, 

the  power  spectral  density  of  n-th  finite  difference,  S . . (co),  of  the 

(Ax) 

function  X is  related  to  the  power  spectral  density  of  X by  the  approximate 
relation 


S 


(«) 


SX(U)  • (UT 


Thus,  if  S-^(co)  is  low  frequency  divergent  in  a power-law  sense, 
there  always  exists  some  finite  n for  which  the  n-th  finite  difference  has 
a convergent  power  spectral  density.  This  fact  is  exploited  by  Barnes 
[1966]  for  the  recognition  and  classification  of  flicker  noise.  Provided  n 
is  sufficiently  large,  the  limit 


Lim  N 
N-* 00 


(2-3) 


always  exists.  As  before,  if  the  measuring  system  has  a high  frequency 

cutoff,  co  , then  one  normally  works  in  the  range  where 
h 


co  t > > 1 . 
h 

The  mean  square  of  the  n-th  finite  difference  is,  in  general,  a function 
of  the  delay  time,  t,  and  the  high  frequency  cutoff,  co^. 

It  is  true  that  if  N (the  number  of  independent  differences  used)  is  not 
large  or  the  order  of  differencing,  n,  is  not  adequately  large,  the  expected 
value  of 

N 

-nZ 

i=l 

for  fixed  N may  differ  from  ^AnX^^  as  defined  in  (2-3).  One  may,  in 
fact,  use  Allan's  [ 1 96 6]  theory  to  treat  these  cases  (see  also  Section  2.  3). 
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2.  6 Ratio  of  Variances  of  Finite  Differences 


In  analogy  to  the  X-function  of  Allan  [1966]  (see  Section  2.  4),  one 
may  compute  the  ratio  of  the  variance  (or  mean  square  if  appropriate) 
of  the  (n  + l)-th  finite  difference  to  that  of  the  n-th  finite  difference.  This 
ratio  is  a function  only  of  the  order  of  differencing  used  and  the  noise 
spectral  type.  By  extending  the  work  of  Barnes  [ 1 966] , figure  5 can  be 
obtained  for  the  ratio  of  variances  of  finite  differences.  As  shown  in 
figure  4,  the  point  at  a = 1 is  actually  degenerate  for  all  a such  that  a ^ 1. 
One  might  consider  that  for  this  situation  ^ is  the  more  fundamental 
quantity  and  fj,  ^ - 2 as  indicated  in  figure  4 regardless  of  a. 


2.  7 Bandwidth  Dependence 

Vessot,  et  al.  [1966]  consider  the  case  where  one  may,  roughly, 
equate 

1 

T *3  


and  consider  <j 2 I00, — , 0,co  ).  Obviously  this  will  exist  only  for  noises 
\ uh  h/ 

where  a < 1 . This  has  been  referred  to  as  the  "marching  bandwidth" 

method.  As  one  may  expect  there  does  exist  noise  spectral  types  where 

the  variation  of  <j  with  or  is  one  of  the  better  means  of  diagnosis. 

h 


2.  8 Summary  of  Relations 

The  use  of  precision  oscillators  in  frequency  standards,  clocks, 
spectrum  analyzers,  and  Doppler  radar  ranging  has  resulted  in  the 
development  of  several  statistical  quantities  of  significance.  By  virtue 
of  the  variable  use  of  oscillators  it  is  of  value  to  relate  the  various 
quantities  and  consider  their  sensitivity  as  a diagnostic  tool. 

Consider  §(t)  to  be  the  fluctuations  in  phase  of  some  oscillator 
relative  to  an  ideal  definition.  The  instantaneous  frequency,  <J>  (t)  = — <}>(t), 
is  assumed  to  have  a power  spectral  density,  S • (oj),  which  may  be 
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kD 


<5(xuv>> 

<=<XI+u^» 


'oi«H 


32 


Exponent, 


expanded  in  a Laurent  series  where  only  four  terms  are  of  any  significance 
in  its  use.  That  is, 


s-  M 

9 


A I co  I + B + C I co  I + Deo  , I co  | < oo 

h 

0,  otherwise. 


(2-4) 


This  specifically  assumes  that  the  lowest  order  noise  process  affecting 
the  frequency  is  a flicker  noise.  In  actual  practice  if  the  oscillator  is 
subjected  to  various  external  perturbations  one  would  expect  a term 
proportional  to  |oo  | “2  to  be  of  significance.  The  present  model,  thus, 
assumes  "ideal"  ambient  conditions. 

It  is  of  value  to  consider  four  different  statistics  of  the  process: 

First,  the  standard  deviation  of  the  average  frequency,  cr2  ^N,  t,  t^,  aj^j, 
as  defined  in  Section  2.  3 is  important.  Second,  when  the  dead  time,  t^, 
is  zero,  the  equations  simplify  significantly  and  it  is  of  value  to  consider 
this  special  case.  The  third  quantity  of  interest  is  the  mean  square  second 
difference  of  the  phase,  D2(t,g)  ) = ((A2^)2). 

When  used  in  Doppler  radar  ranging,  one  compares  the  phase  elapsed 
during  some  interval,  t,  with  that  elapsed  during  a delayed  interval  of 
the  same  duration.  Thus,  the  fourth  quantity  may  be  defined  as 


ip  c(t,  T,  u)  ) - ([cj>(t  + T + t)  - cf>(t+T)  - cj>(t  + t)  + 4>(t)]2  > 


and  will  be  referred  to  as  the  Doppler  range  error.  Table  2 lists  the 
definition  of  the  four  quantities  and  relates  them  to  the  autocovariance 
functions  and  the  power  spectral  densities  for  general  (effective 
stationary)  noise  processes. 

Table  3 shows  general  relations  between  the  four  statistical  quantities 
for  any  noise  process.  For  the  specific  model,  (2-4),  the  terms  arising 
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from  each  term  of  the  Laurent  expansion  are  given  for  each  of  the  four 
statistical  quantities. 

Inspection  of  Table  4 reveals  various  general  rules  for  diagnostics  of 
power-law  spectra: 

(1)  Changing  N in  the  variance  (X-test  of  Section  2.  4)  is  useful  in 
recognizing  white  or  flicker  frequency  modulation  (A  and  B 
terms)  but  not  for  distinguishing  between  white  and  flicker 
phase  modulation  (C  and  D terms). 

(2)  Changing  t for  any  of  the  quantities  is  similar  to  changing  N 
((1)  above)  in  the  variance. 

(3)  Changing  is  the  most  sensitive  method  of  distinguishing  white 
phase  noise  (D-terms)  from  the  others  (A,  B,  and  C-terms). 

(4)  Increasing  from  zero  in  the  variance  has  the  general  result  of 
losing  sensitivity  for  distinguishing  among  any  of  the  terms. 
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G eneral  Relation 
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LAURENT  EXPANSION  OF  SPECTRAL  DENSITY 
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Appendix 


A.  1 The  Random  Walk  Problem 

Consider  a filter  whose  transfer  function  may  be  expressed  as 


gH  = 


(A-l) 


where  oo^  is  the  low  frequency  cutoff.  It  is  easy  to  realize  this  exact 
transfer  function  with  a resistor  and  a capacitor. 

As  in  Section  1.  1,  assume  the  input,  e^,  is  a white  noise  whose 
spectral  density  is 

S.(oo)  = h. 

l 

The  output,  e^,  has  a spectral  density  which  is,  then. 


SQ(U))  = |g(w)  |aS.(u), 


(A-Z) 


or 


s>)  = I ps 

luf  + jwl 


(A-  3) 


From  work  with  analog  computers,  one  knows  that  as  becomes 

smaller,  the  filter  becomes  more  nearly  a perfect  integrator.  Thus  one 

h 

is  tempted  to  consider  the  limit  function  of  the  sequence 


!- 


as  co 


i 


|oo,g  + jtt  fj 

0 to  be  the  true  power  spectral  density  of  a (continuous)  random 
walk.  This  limit  function  is  just  h | co  | 2. 

For  non-zero  oo^,  one  may  obtain  the  autocovariance  function  of  the 
output  defined  by  the  relation  T 


R (t)  = Lim 
o m 

T-*  00 


7T  f e (t)e 
T / o o 


(t  + f)dt. 


(A-4) 
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4 


It  is  often  easier  to  obtain  R (t)  as  the  inverse  Fourier  transform  of 

o 

S (go)  (Wiener-Khinchin  theorem).  For  u> , / 0 no  problem  is  encountered 
o f 

because  eQ(t)  is  a stationary  time  function.  The  result  is  [Campbell  and 
Foster,  1948]  (pair  444) 

Ro(t)=IT"  exP  f UI  I T l)  • (A~5) 

i 

It  is  apparent  now,  however,  that  the  sequence  {Rq(t)}  does  not  have  a 
limit  function  as  go^  -*  0.  While  one  can  consider  the  limit  of  the  spectral 
density  sequence,  the  limit  autocovariance  function  does  not  exist,  and, 
thus,  in  the  limit,  the  Wiener-Khinchin  theorem  does  not  apply.  With 
the  autocovariance  not  defined,  one  has  difficulty  working  in  the  time 
domain. 

There  is,  however,  an  interesting  approach  to  this  problem.  One 
knows  that  a random  walk  has  "stationary  increments"  [Parzen,  1962]. 
Assume  to^  ^ 0 temporarily.  One  may  consider  the  quantity 

Ae  (t)  = e (t+r)  - e (t), 
o o o 


and,  thus. 


(Km)^  2 


ro(°) 


(A—  6) 


where  the  angular  brackets  denote  time  average,  and  use  has  been 
made  of  the  definition  of  the  autocovariance  function,  (A-4).  Substitution 
of  (A-5)  into  (A-6)  yields 


(A-7) 


One  may,  now. 


consider  the  limit  of  (A-7)  as  go 


f 


0.  This  limit  is 


Lim 

go  -*  0 
l 


Zv  h t 


(A-  8) 
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The  usual  factor  of  Zv  is  present  because  S (co)  is  defined  to  be  a 

o 

density  relative  to  angular  Fourier  frequencies.  Equation  (A-8)  does,  in 
fact,  give  the  true  mean  square  value  of  the  increment  in  a random  walk 
or  Brownian  Motion  [Parzen,  1962;  and  Cutler  and  Searle,  1966]. 

Some  insight  may  be  gained  by  considering  the  Taylor  series 
expansion  of  (A-5);  i.  e.  , 


R (t)  = Trh 


(A- 9) 


It  is  apparent  that  the  summation  on  the  right  of  (A- 9)  vanishes 
limit  as  o -*  0.  The  term  - t is  independent  of  co  as  go  „ -*  0. 

apparent  that  if  functionals  of  R (t)  are  considered  whose  form 

N 


F(Ro(T>)£>nRo<Tn>- 


n=  1 


in  the 
It  is 
is 


then  the  limit  as  co  „ 

£ 


0 will  always  be  given  by 


N 


Lim  F(R0(t)>  = -£XlTn 


ir  h 


co  -*  0 
£ 


n=l 


N 


provided  ordy  that£Pn  = 0.  With  this  constraint  satisfied,  one  tnay 


n=l 


consider  the  "effective,  stationary  (E-S)  autocovariance  function"  of 

S(co)  = h | co  | “2  to  be  R (t)  = -irh  |t  | . 

ES 

One  may,  therefore,  treat  a random  walk  as  a stationary  process 
and  use  the  E-S  autocovariance  functions  in  the  time  domain  to  solve  the 
resulting  time  domain  problems.  One  must,  however,  be  certain  that 
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* 


the  constraint,^  ^(3^  = 0,  is  satisfied.  A significant  and  interesting  point 


is  that 


res(t)  " 71  |t 


is  the  (generalized)  Fourier  transform  of  the  generalized  function  | co  | 
[Lighthill,  1962].  It  is  obvious  that  E-S  autocovariance  functions 

do  not  necessarily  satisfy  the  usual  inequality 


R (0)  2 |r(t)|. 
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A.  Z Extension  of  Effective  Stationarity 

Consider  the  power  spectral  density  function 


S(w)  = h | co  + jco 


- a 


(A-l  0) 


The  autocovariance  function,  R(t),  corresponding  to  (A-10)  may  be 
found  using  the  Wiener -Khinchin  theorem  and  pair  569  of  Campbell  and 
Foster  [1948J.  That  is  qj.x 


K (w.  T 
a- 1 


R(t)  = Z-rth 


i 


or  - 1 


(A-l  1) 


where  K ^(X)  is  the  Bessel  function  of  the  second  kind  for  imaginary 
argument.  Two  significant  relations  are 


V°  =tA — [Wx) " Vx) 

u Z sm  1/tt  1 


(A-12) 


and 

“ p+  2k 

Vx)  =£  

k=0VZ^  j(k!)r(^  + k+l) 

From  (A-l  3)  the  term  corresponding  to  k = 0 for  V ~ 


O'-  i 


and 


(A- 13) 


X = ^ M is 


(“i  |t  [)" 


<y  - 1 


Substitution  in  (A-1Z)  and  then  (A-ll)  gives  a term  in  R(t)  which  is 

. ^ . a-l 


res(t)  = 


-IT  yfThlj, 


_/or\  _ /ar+ i\  . fTT(o--i) 

r(?)rv— ) sin 


(A-l  4) 


4Z 


♦ 


If 


a - 1 


is  an  integer,  then  the  term  in  (A- 13)  for 


<x  — i a - i 

k = — - — and  V - — also  contributes  to  the  -independent  term  but 

a - i 

— 2 — being  an  integer,  the  denominator 


2 — v*  „ 2 

with  inverted  sign.  However,  with 


of  (A-14)  also  vanishes.  Always  the  special  cases  of  integer  values  of 
ex  - 1 

— - — (flicker  noise  and  its  integrals)  must  be  treated  separately. 

For  a = 2 in  (A-14)  one  obtains,  in  fact,  R (t)  = - Trh  |t  | in 

ES 

agreement  with  Section  A.  1.  For  a being  an  even  integer,  (A-14)  reduces 
to 


SL  a - i 

2 |T  I 

RES  = <_1)  1111  (a-1)! 


(A- 15) 


which  is  in  exact  agreement  with  the  Fourier  transform  of  the  appropriate 
generalized  functions  [Lighthill,  1962].  Making  use  of  the  duplication 
formula  [Whittaker  and  Watson,  1952]  for  gamma  functions, 


22Z_1r(z)r(z+  i)  =^r(2z), 

(A-14)  may  be  written  in  the  form 


Res(t)  = irh  |t 


a - i 


1 


r (q')cos 


(A-16) 


which  is  also  in  exact  agreement  with  generalized  functions.  Equation 

(A-16)  does  not  apply  for  or  being  an  odd  integer,  however. 

<x  — l 

The  special  cases  of  — - — being  an  integer  may  be  evaluated  as  the 
limit  of  the  sum  of  the  two  important  terms  in  the  expansion  of  (A- 11)  as 

O'  -1 

— - — approaches  the  integer  n.  The  first  term  is  just  that  given  by  (A-14). 
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A.  3 Constraints  on  Effective  Stationarity 

Consider  a function,  <j)(t),  and  a second  function,  ^j(t),  obtained  from 
4> (t ) by  the  relation 

N 

i/l(t)  =^bn4>(t  + CnT)  (A- 1 8) 

i=l  n 

where  the  set  [b  , C } and  t are  constants.  If  one  is  interested  in  cal- 
n n 

culating  the  mean  square  of  )p(t),  this  can  always  be  expressed  in  the 
form 

M 

<(<Mt))2>  =£piiyYiT)  <a-i9> 

i=l 

where  the  set  {(3.,  7.}  are  obtained  from  the  [b  , C } of  (A-18)  and  it  is 
ii  n n 

assumed  that  <f>(t)  is  effectively  stationary. 

It  has  already  been  shown  that,  for 


s , (<*>)  = h |u  | 2, 

9 

the  constraints  on  are  given  simply  by  the  relation 

M 

I>i  = °- 

i=l 


It  is  a simple  matter  to  show  that  the  condition 

N 

Ebn=° 


n=l 


for  the  coefficients  [b  ,C  } is  also  sufficient  to  guarantee  correct  results. 

n n 

As  shown  by  Barnes  [19  66],  for  a ^ 3,  the  constraints  on  [b  , C } 

n n 

become 


= 0 
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and 


N 

Zb  b (C  -C  f = 0. 
nil  n 

n<! 

Intuitively  one  would  expect  these  constraints  to  break  down  at  a 2.  5 
since  a second  difference  does  satisfy  the  constraints  but  fails  to 

converge  for  a 2 5 (see  [Barnes  and  Allan,  196?]).  These  conditions 

, M 

become  t — > 

i=  1 


and 


M 

- 0 
i=l 


for  the  {(3.,"Y.}.  By  an  extension,  one  would  expect  that  the  constraints 

M 

I>i  = 0 

i=l 

M 

"° 

i=l 


M 


Y , = 0 


Epi-i 

i=l 

should  guarantee  effective  stationarity  to  the  region  a < 7.  One  may 
quickly  verify  that  the  third  finite  difference  satisfies  all  of  these  con- 
ditions. Similarly  one  could  add  the  condition 


N 


Vb  b (C  -C  )4  = 0 
L~i  nil  n 

n <1 

for  the  [b  , C }.  These  results  are  compiled  in  Table  6. 
n n r 
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TABLE  6 

COEFFICIENT  CONSTRAINTS  FOR  EFFECTIVE  ST ATIONARIT Y 


Range  for  a 

{b  C 1 Constraints 
n n J 

{Piy^}  Constraints 

S+M  = h |w  | 01 

N 

m =Ebn^(t+CnT) 

n= 1 

M 

V <t»  =X)PiR4»'v  iT) 

i=  1 

a < 1 

None 

None 

(<(>  is  stationary) 

( c|)  is  stationary) 

a < 3 

N 

Vb  =0 

l n 
n=  1 

o 

ii 

• i— < 

sRE 

a < 5 

N 

= 0 
N 

Vb  b (C  -C  )2  = 0 
Z— i n i i n 

n<i 

M 

I>i  = 0 

i=l 

M 

&.v.s  = 0 

i=  1 

a < 7 

N 

= 0 

n=  1 n 
N 

V b b (c  - C f = o 
Zi  n r 1 n 

n<i 

N 

Vb  b (C  - C )4=  0 
Z-/  n 4.  !.  n 

n<i 

M 

&■-" 

M 

I>iV  = ° 
1=  1 

M 

= 0 

iti  1 1 
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